Foldy-Wouthuysen transformation for relativistic particles in 

external fields 



Alexander J. Silenko 

Institute of Nuclear Problems, Belarusian State University, 
Bobruiskaya str., 11, 220080 Minsk, Belarus 
(Dated: February 5, 2008) 

Abstract 

A method of Foldy-Wouthuysen transformation for relativistic spin-1/2 particles in external 
fields is proposed. It permits determination of the Hamilton operator in the Foldy-Wouthuysen 
representation with any accuracy. Interactions between a particle having an anomalous magnetic 
moment and nonstationary electromagnetic and electroweak fields are investigated. 

PACS numbers: 03.65.Pm, ll.lO.Ef, 12.20.Ds 



1 



I. INTRODUCTION 

The Foldy-Wouthuysen (FW) representation [l| occupies special place in the quantum 
theory. This is mainly due to the fact that the Hamiltonian and all operators in this repre- 
sentation are block-diagonal (diagonal in two spinors). For relativistic particles in external 
fields, operators have the same form as in the nonrelativistic quantum theory. Therefore, 
the FW representation in the relativistic quantum theory is similar to the nonrelativistic 
quantum theory. The basic advantages of the FW representation are described in see 
also below). 

The transformation to the FW representation (FW transformation) holds only in the 
one-particle approximation where the radiative corrections are not calculated in a consistent 
way but are phenomenologically taken into account by including extra terms in the Dirac 
equation (see ^). One-particle description is feasible even for ultrarelativistic particles if 
the external field is so weak that the probability of pair production or bremsstrahlung losses 
can be neglected for a given interaction energy of a particle. The range of applicability of 
this description is fairly wide and includes, in particular, the relativistic particle scattering 
and the interaction of relativistic particles with matter and external fields. 

In the nonrelativistic case, there exist a lot of good methods of FW transformation with 
taking into account relativistic corrections However, they are not useful for 

relativistic particles. The known methods of solving this problem Q, U,^^ either lead 
to cumbersome calculations or the field of their use is limited by the first approximation in 
field parameters. None of these methods permits exact FW transformation for the particular 
cases described in . Therefore, the FW representation does not take the right stand 

in the relativistic quantum theory. The Dirac and Melosh representations are mostly 
used. 

FW transformation can also be performed for particles with spin s > 1/2 j?! 

In the present work, a method of FW transformation for relativistic particles in external 
fields is proposed. This method permits obtaining a Hamiltonian of any accuracy by succes- 
sive approximations, as a power series in the external field potentials and their derivatives. 
In some cases, this method permits performing exact FW transformation. 

The relativistic system of units h = c = 1 is used. 
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II. GENERAL PROPERTIES OF THE FOLDY-WOUTHUYSEN REPRESENTA- 
TION 



The basic advantages of the FW representation are due to its specific properties. 

The relations between the operators in the FW representation are similar to those be- 
tween the respective classical quantities. In this representation, the operators have the same 
form as in the nonrelativistic quantum theory. Only the FW representation possesses these 
properties considerably simplifying the transition to the semiclassical description. The FW 
representation provides the best possibility of obtaining a meaningful classical limit of the 
relativistic quantum theory. 

For example, the Hamiltonian for a free particle fully agrees with that of classical physics: 



Hfw = i5\/m^^p', p = -^V, (1) 

in contrast with the Hamiltonian in the Dirac representation ^, 15|. The position operator 
in the Dirac representation is the radius- vector, r [15]. It corresponds to the mean position 
operator for the free particle in the FW representation 



i3a i3(cx ■ p)p + lYl X p\p 
rD = r + — — — ■ ^ , p 



\P\ 



2e 2e(e + m)p 

Here and below the following designations for the matrices are used: 



7 



cr 
-cr 



OL = 






1 
-1 

n = /3s 



where 0, 1, —1 mean the corresponding 2x2 matrices and cr is the Pauli matrix. 

In the FW representation, the problem of zitterbewegung'^ motion never arises 
The operators I = r x p and S/2 define the angular momentum and the spin for the free 
particle, respectively. In this representation, unlike the Dirac one, each of them is a constant 
of motion (see ^). The corresponding operators conserving in the Dirac representation are 

= r-D X p, 
_ S tf3[cx xp] [p X [S X p]] 
2 2 2e 2e(e + m) 
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The total angular momentum operator, ji, is a constant of motion in both representations, 
because 

Jd = Id + — =1 + ^=J- 
The FW representation is very convenient for describing the particle polarization. In 



this representation, polarization operators have simp 
dimensional polarization operator equals the matrix U 



this operator depends on the particle momentum 

o = = n - 7 



3, 



e forms. For example, the three- 
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7[ . In the Dirac representation. 



e e(e + m) 

i^'or particles interacting with external fields, it also depends on the external field parameters 



Thus, in the Dirac representation all operators corresponding to the basic classical quan- 
tities are defined by cumbersome expressions. These operators should also depend on the 
external field parameters for particles interacting with external fields. 

The FW representation helps one to prove that the particle position can be measured 
up to its Compton wavelength Q]. However, this property is valid only for a particle 
not strongly interacting with external fields if the one-particle approximation is attainable. 
Otherwise, the effect of pair production prevents the use of both the Dirac equation (even 
with some corrections) and the "traditional" Hamilton approach. Obviously, in this case 
the FW transformation cannot be used either. 

The FW transformation possesses another important property. The relativistic wave 
equations and all operators are block-diagonal (diagonal in two spinors). This property 
permits separating positive and negative energy states Q|. Of course, extraction of even 
parts of operators becomes unnecessar y. 

The detailed analysis performed in [1^ shows that the wave functions in both the Dirac 
and FW representations are equal to each other only approximately, and they do not co- 
incide. In these papers, the nonrelativistic case was considered and relativistic corrections 
was taken into account. 

An analogous conclusion follows from the results obtained in In this work, a more 



general situation has been investigated for the relativistic particle not strongly interacting 
with an electromagnetic field. It has been found that the upper spinors in the Dirac and 
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FW representations are approximately proportional to each other, but this property is not 
exact. 

Thus, the preferable employment of the FW representation is evident, although the rel- 
ativistic wave equations are more complicated in this representation. 



III. METHODS OF THE FOLDY-WOUTHUYSEN TRANSFORMATION 

In the classical work by Foldy and Wouthuysen P], two different transformations, for 
free relativistic particles and for nonrelativistic particles in electromagnetic fields have been 
carried out. In the general case, transformation to a new representation described by the 
wave function \E'' is performed with the unitary operator U: 



where = | | is the wave function (bispinor) in the Dirac representation. As 
X ^ 

at ot 

the following transformation can be carried out: 



I 



dt / ' \ dt 



Hence, the Hamilton operator in the new representation takes the form jl,|l9|: 



or 



n' = unu-^ - tu^—, (2) 

ot 



u' = u \ n-i—\u 



i- 



dtj dt' 
There is an error in this transformation in |20]. 

The Hamiltonian can be split into operators commuting and noncommuting with the 
operator /?: 

n = (3m + S + 0, f3S = S(3, (30 = -0(3. (3) 
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The Hamiltonian Ti is Hermitian. We assume that both operators £ and O are also 
Hermitian. 

For free Dirac particles £ = 0, O = ex. ■ p, and the operator S has the form 



S = -ipcfp9{p), 



(4) 



where ^ is a function of the momentum operator. If we choose 

^(p) = — arctan f — ) , 
2p \m J 

the transformed Hamiltonian H' contains no odd operators |l, 20| and we obtain Eq. (1) 



For nonrelativistic particles in an electromagnetic field, the FW transformation can be 
performed with the operator 



The transformed Hamiltonian can be written in the form 

n' = n + z[s, n] + '^[s, [s, n]] + '^[s, [s, [s, n]]] + 



(5) 



(6) 



where 



means a commutator. As a result of this transformation, we find 



n' = f3e + S' + O', (38' = S'(3, (30' = -0'(3, 



(7) 



where the odd operator O' is now 0(l/m). This procedure can be repeated to obtain the 
required accuracy. Another form of the nonrelativistic FW transformation was given by 
Ericksen (see also 21|). 



There are also other methods for obtaining the block-diagonal form of the Hamiltonian 

n 

or the Lagrangian. The so-called elimination method of Pauli [22] permits excluding the 
lower spinor from relativistic wave equations. As a result, the wave function of the final 
Pauli equation is the upper Dirac spinor, 0. This means that the upper Dirac spinor is also 



an eigenfunction of the transformed Hami 



Pauli method was analyzed in detail in 



tonian. However, this property is not exact. The 



21|. It was shown that this method gives the right 



first approximation. Nevertheless, relativistic corrections of higher orders are incorrect. It 



is quite natural because direct Pauli's reduction leads to a neglect of the contribution of the 
lower spinor The relation between the exact wave function in the FW representation 
and the upper Dirac spinor has been found in in the relativistic case. 

A more exact variant of the elimination method had been proposed earlier by Berestetskii 



and Landau j23| (see also |2^). They showed that it was possible to find a nonunitary 
operator V for which 

ij = V(j) (8) 

is a two-component wave function with a correct norm. An appropriate form of the operator 
V can be obtained from the condition 



The relation between the Dirac spinors can be expressed in the general form: x = Q4^- 
Therefore, 

V^V = 1 + Q^Q. 

If we additionally assume that the operator V is Hermitian, then both this operator and 



23 



2i 



25|. 



the Hamiltonian can be found by successive approximations 

Of course, the elimination method is much simpler. However, it is mostly intuitive. Its 
validity is proved only by the coincidence of the results obtained by the FW and Akhiezer- 
Berestetskii-Landau methods |2l| . 

Another method of diagonalization of relativistic wave equations was proposed by Korner 
and Thompson ^. In this work, the Lagrangian approach was used. The Korner-Thompson 
method is similar to the FW method. It also includes a successive decrease in the maximum 
order of odd terms. The results obtained by the FW and Korner-Thompson methods agree 
(see Q). 

Thus, several nonrelativistic transformation methods give the same results. However, the 
FW transformation method has been justified in the best way. 
In several cases, FW transformation can be performed exactly 

urn 

. Exact FW trans- 
formation has also been performed for a wide class of external fields in j^]. In this work, 
involutive symmetries of relativistic wave equations have been used. However, the trans- 
formed Hamiltonians contain "nontraditional" space reflection operators. The reduction of 
Hamiltonians to the "traditional" form is a difficult problem. It has not been investigated 
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in 13]. However, this reduction is necessary to do for solving many problems (e.g., finding 
particle and spin motion equations). 

Generally, FW transformation for relativistic particles in external fields is complicated. 
The transformation methods explained in ^ require cumbersome calculations. A variant 
of the elimination method useful for relativistic particles has been developed in On 



eliminating the lower spinor 
upper spinor takes the form 



rom the relativistic wave equations, the final equation for the 



d(j) d ^ , 



where F is the operator function. Further calculations are analogous to those in the 
Akhiezer-Berestetskii-Landau method. A new wave function with a correct norm, ip, ex- 
pressed by Eq. (8) is introduced. Substituting it for into Eq. (9), one can find the 
Hamilton operator for the relativistic particle. 

The relativistic wave equation for an upper spinor similar to Eq. (9) is found by the 



Lagrangian approach 
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However, it is difficult to find a second approximation by using the relativistic variant of 
the elimination method proposed in It is easier to determine relativistic corrections of 

n 

higher orders |12l |. 

The right two-component wave function in the FW representation, ip, does not coincide 
with the upper Dirac spinor, ^|. This conclusion is in agreement with the results obtained 
in Q- 

There are other difficult problems. The diagonalization of relativistic wave equations 
needs carefulness, especially in the time-dependent case. As mentioned in Q, |2^, in the 
latter case Ti' is not equivalent to Ti. since these operators have different matrix elements. 
Rather, UHU'^ is. There is a danger that one can arrive at a block-diagonal representation 
differing from the FW one even in the time-independent case. For example, the transforma- 



«on perfo.n,ed . ,tM. . the Me,o.h .a.fo..atio„ .deed M) lead= to a bIoc.-dia.onal 

Hamiltonian that differs from the Hamiltonian in the FW representation [2^. Therefore, 
the application of noncanonical transformation methods is restricted by the necessity of 
verifying the results by comparing them with the corresponding results obtained by the 
canonical transformation method in some particular cases. Of course, other transformation 
methods may be simpler or less cumbersome. Nevertheless, the FW method is safer and 
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substantiated very well. 

In the present work, a relativistic extension to the FW method is proposed. 



IV. EXACT FOLDY-WOUTHUYSEN TRANSFORMATION 

Consider some cases of the exact FW transformation. 

In Eq. (3), the operators (3 and commute {f30^ = -Of30 = 0^f3). Therefore, the 
operator is even. 

The operator S can be defined by an expression similar to Eq. (4): 

S = -^^e, (10) 

where C and 9 are the functions of and the operator C satisfies the following conditions: 

C' = C»^ [P,C] = 0. (11) 

It follows from conditions (11) that the operator C is also even. 

It is possible to use the following formal definition of this operator: 

C = VO^. (12) 

Relations (11), (12) define the square root of matrix operators. To unambiguously define 
the square root, these relations should be complemented by the condition that the square 
root of the unit matrix I is equal to the unit matrix. This definition of the square root 



coincides with those of 



201 . |26[ . For example, for free particles 



= ap, 0^=lp\ C = 1^ = I\p\. 

Further, the symbol of the unit matrix I will be omitted. 
Since 

= -(30(30, C = ^-(30(30, f{0^) = f{C^), 

the operators (30^ C, O'^, and 9 commute with each other. The operator 9 is the angle of 
rotation of the basic vector set in the spinor space. 

The unitary transformation operator, U, can be written in the form 

BO 

U = cos9 + ^sm9. (13) 



An FW transformation is exact if the external field is stationary and the operators 8 and 
O commute: 

[£,O] = 0. (14) 

In this particular case, 

Condition (14) is a sufficient but not necessary condition of the exact transformation. 
The Hamilton operator in the new representation takes the form 

W = (^cos9 + f^sm9^ n l^cos^- ^sin^ 

= ipm + O) (^cos9 - ^sm9^ + S = {Pm + O) (^cos 2^ - ^ sin 2^^ +S 

cos 2^ -— sin 20 j +S. 

The Hamiltonian Ti' is even if the odd term (proportional to O) vanishes. This takes 
place if 

C 

tan 20 = —. (15) 

m 

This equation has two solutions, 9i and 02, differing in 7r/2. Since 

2tan0 ^ tan20 

tan20= tan0 



l-tan20' 1± Vl + tan2 20' 

they are defined by the relations 

C „ C 



tan0i = , tan02 = , e = ^m? + = ^m? + O^. (16) 

e + m e — m 

Thus, there are two unitary transformations of the operator Ti to an even form. They 
are characterized by the angles 0i and 02, where the angle 0i corresponds to the FW trans- 
formation. 

As a result of both transformations, one of the spinors (lower for 0i and upper for 02) 
becomes zero as for free particles. 

Note that the transformation under consideration is also similar to the Melosh transfor- 
mation 

If condition (14) is satisfied, then the Hamilton operator in the FW representation is 
defined exactly: 

n' = f3e + £. (17) 
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Unlike j9|, |10|, |ll|, Eq. (17) contains exact expressions for the Hamiltonian derived in 



as particular cases. 

The transformation operator U can be written in nonexponential form. After the calcu- 
lation of sin^^i and cos^i with formulae (16), 



sin6'i = — =, cos6'i = \ , 

v/2e(e + m) ' V 2e ' 

we obtain the following expression: 

u- = i±^±m. (18) 

v/2e(e + m) 

where U+ = U, = U-\ This expression agrees with the well-known formula for free 
particles Since [jSOY = 0(3 = —(30, the operator U is unitary. Simultaneous change 
of signs of sin6'i and cos^i does not affect the final result because the wave functions \1', 
are determined up to a sign. Direct calculation of the Hamilton operator in the FW repre- 
sentation also leads to Eq. (17) in accordance with formulae (2), (3), (14), (18). 

Another class of Hamiltonians permitting exact FW transformation has been investigated 

in ha. 



V. EXACT TRANSFORMATION FOR PARTICLES IN ELECTROWEAK 
FIELDS 

Let us consider the interaction of a relativistic spin-1/2 particle, possessing an anomalous 
magnetic moment (AMM), with stationary electromagnetic and electroweak fields. The 
Hamiltonian of the electromagnetic interaction is defined by the Dirac-Pauli equation [s^ : 

HDP = Ci-7T + (3m + e^ + l^'{-Il- H + i^-E), n = p - eA, (19) 

where /i' is AMM, $, A and E, H are the potentials and the strengths of an electromagnetic 



field. This equation is derived in the one-partic 



e approximation and is useful when an 



nl). 



electromagnetic field is not extremely strong (see 

The weak interaction Hamiltonian should be added to the Hamiltonian of Eq. (19). The 
weak interaction does not conserve the spatial parity. For the interaction transferred by neu- 
tral currents, the standard model gives the following expression for the parity-nonconserving 
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weak interaction Hamiltonian in the approximation of a small transferred momentum 

V-PNc = + C2CX ■ a') n(r), 7^ = |^ j , (20) 

where G is the Fermi constant, cr' is the Pauli matrix for matter particles, and n{r) is the 
density of matter particles. For the interactions with nuclei, n{r) characterizes the density 
of nucleons of a certain kind, and cr' should be replaced by the nucleus spin. Formulae 
(19), (20) do not change if the external fields are nonstationary. The matter particles are 
considered to be at rest. 

The coefficients Ci, C2 are different for different pairs of interacting particles. The Hamil- 
tonians corresponding to the interactions with different matter particles should be summed- 
up. The signs in formula (20) depend on the definition of the coefficients Ci, C2 and matrix 
7^. The total Hamiltonian equals 

'H = Hdp + Ti-PNC- (21) 
In this case, in formulae (3),(16)-(18) we have 



£ = - IJ,'U- H, = a-7r + tfi'j ■ E - {Ci-f^ + ■ a') n{r) 

V2 



P 



G 

7 ■ TT + ifi'a -E- — {GiP-f^ + G2I ■ (t') n{r) 
v2 



(22) 



Let us consider some particular cases where Hamiltonian (21) satisfies condition (14). 
For these cases the FW transformation is exact. The general case will be analyzed in the 
next section. 

The exact Hamiltonian in the FW representation is given by Eq. (17), where £ is defined 
by Eq. (22), and 

1/2 



G 

7 ■ TT + i/i'a ■ E-—= (Ci/97^ + G2I ■ cr') n{r) 
V 2 



+ TT^ + /3/i' (S ■ [-KxE] - S ■ [Ex'n:]-V-E) + ^'^E^ -eH ■ H 
+-^(^Ci{S-7r,n(r)}+-C2{o-'-7r,n(r)}+ + C2[Sxo-']-Vn(r) 

-2|3^i'G2[T,xa']■En{r)^ + '^n\r) [GI + SC^ - 2C2(Ci + ^2)2 ■ cr'] 



(23) 
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Hence, the operator is expressed by formulae (18), (22), (23). Although Eq. (23) is 
formally exact, the small terms proportional to C^,CiC2,C2 are wittingly negligible in the 
approximation of a small transferred momentum. 

Formulae (17), (22), (23) describe the exact Hamilton operator in the FW representation 
in the following particular cases: 

a) in the presence of only weak interaction ($ = 0, A = 0, -E = 0, -ff = 0); 

b) for Dirac particles (/i' = 0) in magnetic and weak fields ($ = 0, £^ = 0); 

c) for uncharged particles with AMM in electric and weak fields (e = 0, A = 0, = 0); 

d) for particles with AMM moving in the plane orthogonal to a static uniform magnetic 
field ($ = 0,£; = 0,P, = 0,Ci = C2 = 0); 

e) for uncharged particles with AMM moving in the plane orthogonal to a static 
uniform magnetic field. A static electric field (possibly nonuniform) is also orthogonal to 
the magnetic field (e = 0, ^;±/f, P^ = 0, Ci = C2 = 0). 

In two cases (d) and e)), H = Hsz, and in the case e), the electric field strength does 
not depend on z. Otherwise, rot^^^^O and the magnetic field is not constant {dH /dt^O). 
Therefore, in these cases the operator pz = —i{d/dz) commutes with the Hamilton operator 
and has eigenvalues = const. Consequently, the consideration of the particular case Pz = 
is quite reasonable. All these cases satisfy condition (14). 

Formulae (17), (22), (23) agree with all exact expressions of the operator Tt' derived for 
uncharged particles with AMM in an electrostatic field, Dirac particles in a static magnetic 
field, and particles with AMM moving in the plane orthogonal to a static uniform magnetic 
field in p, ID, ls| • The weak interaction is not considered in these works. 

VI. GENERAL CASE 

In the general case, relativistic particles interact with external fields. We suggest to 
perform the FW transformation in two stages. First, a transformation similar to the FW 
transformation for free particles, is performed for particles in external fields. Second, a 
transformation similar to the FW transformation for nonrelativistic particles is carried out. 

We assume that the external fields are not extremely strong and the transformed Hamil- 
tonian can be expressed as a power series in the field potentials and their derivatives. The 
external fields can be nonstationary. 
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In the general case, formula (17) is not exact because £ depends on the coordinates 
and contains Dirac matrices. We should calculate the commutator of the operators U and 
£-i{d/dt): 



In this case, it is necessary to compute some commutators containing inverse operators 
and square roots of the operators. These commutators can be calculated using the following 
exact formulae which are valid for arbitrary operators A and B 3^ : 

[A-\B] = A-^[B,A]A-'\ (24) 

[A 5] = \{A-\ [A\B]}^ - \ [[A [A,B]IA-'] , (25) 

where A"^ = 1/A and {...,...}+ stands for the anticommutator. If A is the square root of 
the operators and the commutator of the operators is small compared to their product, i.e., 

\[A,B]\ « \ABl 

formulae (24), (25) allow us to obtain the quantity [A, B] with any accuracy by the method of 
successive approximations (see js^). As a rule, this condition is satisfied since it is equivalent 
to the inequality 

I « Ic, (26) 

whereii^ is the total energy including the rest energy and Ic is the characteristic size of the 
nonuniformity region of the external field. For the nonrelativistic particle, the quantity he/ E 
is equal to the Compton wavelength. 

First, it is necessary to perform a unitary transformation with operator (18). After this 
operation, the Hamiltonian Ti' still contains odd terms proportional to the derivatives of the 
potentials. Let us write the operator Ti' as 

n' = pe + 8' + 0', (3S' = S'f3, (30' = -0'(3, (27) 

where 



e = + 0\ 

c>_.d ^ e + m / J\ e + m PO / J\ PO 

dt ^2e(e + m) V dtj ^2e(e + m) v^2e(e + m) V dt J ^2e(e + m) ' 
^, ^ /30 _ e + m e + m _ \ _J30__ 

v/2e(e + m) V ^ 9t J ^2e(e + m) ^2e{e + m) \ ^ dt ) ^2e(e + m) ' 

(28) 
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Since 



ABA = -{{A\B}+-[A, [A,B]]) 



relation (28) for the operator £' takes the form 



8' = 8- 



1 

+ 4 



e + m 



v/2e(e + m) 



e + m 



(30 



8-1 



(29) 



^2e(e + m)' [v^2e(e + m)^ 

The odd terms are small compared to both e and the initial Hamiltonian Ti. This cir- 
cumstance allows us to apply the usual scheme of the nonrelativistic FW transformation 

Second, the transformation should be performed with the following operator: 



U' = exp{iS'), S' 



(30) 



The further calculations are similar to those given in |2fl]. The particle mass should be 
replaced by the operator e noncommuting with the operators 8\ O'. If only major corrections 
are taken into account, then the transformed Hamiltonian equals 



n" 



(31) 



This is the Hamiltonian in the FW representation. 

To obtain the desired accuracy, the calculation procedure with the transformation oper- 
ator (30) (S" is replaced by S" , S'" etc.) should be repeated multiply. 

Let us calculate the Hamiltonian in the FW representation for the relativistic particle 
with AMM interacting with a nonstationary electroweak field. The Hamiltonian in the Dirac 
representation is defined by formulae (19)-(21). The transformed Hamiltonian is defined by 
Eq. (30), where the operator O' contains the field strengths and does not contain the field 
potentials. Let us deduce the Hamiltonian to within first-order terms in the field strengths 
and their first derivatives and second-order terms in the field potentials. The terms of the 
second order and higher in the field strengths and their derivatives and the first-order terms 
containing derivatives of the second order and higher of the field strengths will be omitted. 

Since we neglect the second-order quantities in 0\ the operator O' does not make any 
contribution to the Hamiltonian Ti" at the second stage of transformation defined by for- 
mula (31). As a result, we obtain the following equation for the Hamiltonian in the FW 
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representation: 



n" = f3e + S\ 

£' = e^ + -{—^ — -A1:-\'kxE]-1:-\Extt]-V-E 
+ — <! 7^ ^,7r-V(7r-i; + i;-7r) 1^ -/i'H-if 



32 e4(e + m)2 
1 



(32) 



4 |^e(e + m) 



i3'-7r)(S-7r) + (5]-7r)(7r-if) + 2n{n-j + j-7r) 



1 Si/ 

where j = V x H — is the external current density, and e is determined by formulae 

An ot 

(23), (27). It is important that the operators e,S' are found at the first stage, i.e., at the 
transformation with operator (18). 
In the weak field approximation. 



' ' 1 



4 1 e'' 



and 



4 [ \e + ^ 
fxom ( 2e'^ + 2e'm + 



'E-lExn] - V-E 



+ 



16 
1 



4 1 e'fe' + m) ' 
where 



n.V{n.E + E.n)y-\{{^^ + ^'),U.H}^ 

G (1 

^471 



e'4(e' + m)2 

if-7r)(n-7r) + (n-7r)(7r-i3") + 2n{7r-j + j-TZ 



W 



(33) 



e' = + 7r2, = Ci{S ■ tt, n(r)}+ - Czjo-' ■ tt, n(r)}+ + Cap x ct'] ■ Vn(r), (34) 

and yUo = e/(2m) is the Dirac magnetic moment. 

Ill ], formula (33) includes, as particular cases, the exact expressions 

Formulae (32)-(34) 



Unlike works 



m 



for the Hamiltonian in the FW representation o 



also agree with the results obtained 



in J 



Dtained in 



method of FW transformation used in 



G 
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mm- 

Detailed analysis shows that the 



lll | does not allow one to take into consideration 



the terms proportional to the double commutators of e with e$ and fi'II H (e[e, [e, $]] anc. 
/i'[e, [e, n-i3"]]). The terms proportional to E and H in the Hamiltonian obtained in 11 

and Eq. (33). However, only the Dirac particles were considered 



coincide with those in 
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m lai, the derivatives of the field strengths were neglected 
Hamiltonian with relativistic corrections was found in 



in 0, 



and the nonrelativistic 



VII. PARTICLE AND SPIN MOTION EQUATIONS 



Thus, the FW representation is very convenient for describing the particle and spin 
motion owing to the simple forms of operators. In order to derive corresponding quantum 
equations, it is necessary to compute the commutators of the Hamiltonian with the same 
operators as in the nonrelativistic theory. The kinetic momentum operator of particles in 
an electromagnetic field equals tt = p — eA [35]. The equation of particle motion in the 
electromagnetic field is defined in terms of the commutator of the Hamiltonian with this 
operator: 

dt L ' J dt 

To determine the quantum equation of particle motion, we take into account Eq. (33) 
for the Hamiltonian 



4 I \ e' + m I e' 



V(S ■ [E X tt]) - V(S ■ [tt X E]) + AE 



' -^,{7r-V)Vi'K-E + E-77) 
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4 I e'(e' + m) 



2 



/i' ,v(n-if) 



n-7r)V(7r-ii") + V(/f-7r) (n-7r) + 27rV(7r-j + j-tt 



(35) 



The equation of spin motion is defined by the formula 

For particles in a nonstationary electroweak field, it takes the form 
/iom 



dn 

'dt 



e' + m 

ii{ L_ 

2 le'(e' + m 



+ /i') i, [n X X tt]]}^ + |(/^ + , [S X if]} 
, ( [S X 7r](7r ■H) + {H ■ 7r)[S x tt] 



{c, {[S X 7r],n(r)}^ + [S x [a' x Vn{r)]] 




(36) 
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The corresponding equation for stationary electroweak fields was derived in 3^ . 

The transition to the semiclassical description is also simple. For free particles, the 
lower spinor is equal to zero in the FW representation. For particles in external fields, the 
maximum ratio of the lower and upper spinors is of the first order of Wmt/E, where Wint is 
the energy of the particle interaction with external fields. Thus, we obtain ~ 
(Wint/E)'^. Therefore, the contribution of the lower spinor is negligible and the transition to 
the semiclassical equations is performed by averaging the operators in the equations for the 
upper spinor. It is usually possible to neglect the commutators between the coordinate and 
kinetic momentum operators and between different components of the kinetic momentum 
operator (see 3^). As a result, the operators cr, cr' and tt should be substituted by the 
corresponding classical quantities: the average spin, ^ for matter particles), and the 
kinetic momentum. For the latter quantity we retain the designation tt. The semiclassical 
equations of particle and spin motion are: 



cItt 

'dt 



2V(C ■ [tt X E]) - AE 



fiom 2e + 2e'm + m 



e'4(e' + m)2 



(tt ■ V)V(7r ■ + 



e' 



(37) 



e'(e' + m) 



(|-7r)V(fl"-7r) + 27rV(j-7r) 



^ ^ 2 /^o^ 
dt \e' + m 



2/i' 



+ /i')i[|x[£;x7r]] + 2(^ + /.') [, 



[H ■77)1$, X tt] 



G ( \ 
2C^\i X 7r]n(r) + x x Vn(r)]] . 



e'(e' + m)' v^e' 

Equation (37) shows that the particle motion depends on the spin orientation. The 
corresponding term determines the Stern-Gerlach force. 

It is not as convenient to use the Dirac representation to derive quantum equations 
of particle and spin motion in a similar manner. In this case, it is necessary to extract 
the polarization operator, O, from the obtained equations. This problem is rather difficult 
because the operator O in the Dirac representation is defined by the cumbersome expression 







given m 

For particles in external fields, the FW transformation also changes the form of the kinetic 
momentum operator. In particular, the equation of spin motion in the Dirac representation 
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depends on the operator ttd = U'^irU just as the corresponding equation in the FW 
representation depends on the operator tt. However, these two equations differ in their 
functional dependence on tt. The use of the FW representation protects from both an error 
in derived equations of particle and spin motion and an incorrect interpretation of these 
equations. 

Another method of transition to the semiclassical description is based on the trajectory- 
coherent solution of the Dirac equation [^. 



VIII. DISCUSSION 



As mentioned above, the proposed method permits obtaining a transformed Hamiltonian 
to within first-order terms in the field parameters after the first transformation. In this case, 
an the othe. ca.o„,ca. .nethods need seve.a. t.a,.fon„at,o,. Uflfly Qy There- 

fore, the method described above can be successfully used even for solving nonrelativistic 
problems. For this purpose, the transformation operator, f/, can be expanded in a series 
of 1/m. Of course, such an expansion is helpful only in the case of transformation of the 
operator £ — i{d/dt). The transformation of other operators leads to the appearance of the 
term /3e in Eq. (31). 

Consider the classical example of the FW transformation for a nonrelativistic Dirac par- 
ticle in an electromagnetic field. We calculate the Hamiltonian to within terms of orders of 
{p/mY and p'^W/m^, where W means e$, eA (see jlil^]). In this approximation, the first 
double commutator in Eq. (29) is negligible and the second one is equal to 



[30 



1 ■ TT, 



7- TT, ^ 



2e(e + m 
d 



(30 



m] 



-4 



dt 



[1:-{ttx E)-1:-{Extt)-V ■ E] 



As 



(39) 



e = w? + (a ■ 7r)2 = \/ rn^ + tt^ — eS ■ // = m H 



77 



2m 8m' 

the transformed Hamiltonian is expressed by the well-known formula [Jj, |2 



Si/, 

2m 



n" = (3[m 



2 4 > 



2m 8m" 



6$ 



2m 



U H 



8m2 



[Y. ■ {-K X E) - T. ■ {E X tt) - V ■ E] . 



(40) 
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Thus, the proposed method permits obtaining this formula after the computation of only 
one double commutator. All the other canonical methods require cumbersome calculations 



12 
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2ll |. For example, the classical method of Foldy and Wouthuysen require three 



successive transformations and a calculation of numerous commutators. The noncanonical 
methods (the Pauli's elimination method and others) 



11 



21 



23 



2J, |25| permit 



deriving Hamiltonian (40) in an easier way. Nevertheless, the proposed method is very 
simple even compared to them. Moreover, it gives an opportunity to find a transformed 
Hamiltonian with any accuracy even for relativistic particles in external fields. 

IX. SUMMARY 



In this work, a method of FW transformation for relativistic particles in external fields 



is proposed. This method is sini 
as in known particular cases 



£ and reliable. It performs the exact FW transformation 
as in others. This property distinguishes the proposed 

The method 



method from the other methods developed for relativistic particles jol. Il^ 111 | . 
is based on the well-known elaborations 3, 3- First, a transformation similar to the FW 



transformation for free particles is performed for particles in external fields. Second, a 
transformation similar to the FW transformation for nonrelativistic particles is carried out. 
In the general case, the FW transformation is approximate. As an example, the Hamilton 
operator in the FW representation for relativistic particles with AMM interacting with 
nonstationary electroweak fields is found to within second derivatives of potentials. 
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